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In 2001, the discovery of circumstellar water vapour around the ageing carbon star 
IRC +10216 was announced^. This detection challenged the current understanding of 
chemistry in old stars, since water vapour was predicted to be absent in carbon-rich 
stard^. Several explanations for the occurrence of water vapour were postulated, in- 
cluding the vaporization of icy bodies (comets or dwarf planets) in orbit around the 
staiS, grain surface reactiond^, and photochemistry in the outer circumstellar envelop^. 
However, the only water line detected so far from one carbon-rich evolved star can not 
discriminate, by itself, between the different mechanisms proposed. Here we report on 
the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared 
and sub-millimetre spectrum of IRC +10216, including some high-excitation lines with 
energies corresponding to ~1000 K. The emission of these high-excitation water lines 
can only be explained if water vapour is present in the warm inner region of the enve- 
lope. A plausible explanation for the formation of warm water vapour appears to be 
the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. 
This mechanism triggers also the formation of other molecules such as ammonia, whose 
observed abundanced^ are much higher than hitherto predicted^. 

For stars less massive than about nine times the mass of the sun, the last major nuclear burning 



phase is as an Asymptotic Giant Branch (AGB) stai^^l. A natural chemical division is created between 
C (carbon-rich AGB) stars (with C/0 ratio >1, and hence a surplus of carbon to drive an organic 
chemistry in the envelope, M-type AGB stars (with C/0 ratio <1, yielding the formation of oxygen- 
bearing molecules) and S-type AGB stars (with C/0 ratio ~1). The detection of the ground-state 
transition of ortho- water, ortho-H20(li^o ~ lo,i)i in the envelope around the C-star IRC +10216 came 
as a surpris^, as in thermodynamic equilibrium (TE) chemistry no oxygen-rich molecules (except CO) 
are expected in a carbon-rich environment. The vaporization of a collection of icy bodies (comets or 
dwarf planets) in orbit around the star was invoked to explain the presence of water vapour in this 
carbon-rich environment. It was predicted that water should be released in the intermediate envelope 
at radii larger than (1 — 5) x lO^^crrP. Later on, two other distinct mechanisms were considered as 
possible sources of water vapour observed in IRC +10216, each one making a specific prediction for 
the H2O spatial distribution in the envelope: grain-surface reactions such as Fischer- Tropsch catalysis 
on the surface of small grain^, which would imply that water reaches its maximum abundance at a 
distance around 2 x lO^^cm, and formation in the outer envelope through the radiative association of 
atomic oxygen and molecular hydrogerP. It has also been suggested that water could be formed in 
the warm and dense inner envelop^, although no specific formation mechanism has been proposed 
for such an origin. 

On 12 and 19 November 2009, IRC +10216 was observed with the SPIREP and PAC^ spec- 
trometers on board the Herschel Space Observatorj^Sl. Spectroscopic observations were obtained be- 
tween 55 and 670 fim, at spectral resolving powers between 300 and 450CP3. Currently, many different 
molecules and their isotopologues have been identified: ^^CO, ^^CO, C^^O, H^^qjv^^ H^^qjnj^ ^jj^^ 
SiS, SiO, CS, C^^S, 13CS, C3, C2H, HCl, H^^Cl, ortho-HaO, and para-HsO. The detection of the low- 
excitation ortho-H20(2i^2 — lo,i) transition at 179.5 /im was anticipated, since the energy difference 
between the li^o and 2i^2 level is only 53 K, but the discovery of high-excitation ortho-H20 lines with 
upper level energies around 1000 K came as a surprise (see Fig. 1 and the Supplementary Figs. 1-5). 
These high-excitation ortho-H20 lines provide a strong diagnostic to understand the origin of water 



in carbon-rich envelopes. Moreover, for the first time, para-H20 lines have been detected from the 
envelope of a carbon-rich AGB star (see Fig. 1 and the Supplementary Figs. 1-5). 

The kinetic temperature in the envelope is typically around 2000 K in the dense environment 
just above the stellar photosphere and decreases to ~10K in the tenuous outer envelope (see Fig. 2). 
The presence of high-excitation ortho-H20 lines can only be explained if water is present in the warm 
inner region of the envelope, at radial distances closer than ~15-R* (or 7.5 x 10^^ cm for a stellar 
radius, i?^., of 5.1 x 10^^ cm). This immediately excludes the three mechanisms that only place water 
in the intermediate or outer regions of the envelope as the source of the water origin (see Fig. 2). 
The mechanism involving the radiative association of O and H2 in the outer envelop^ can also be 
ruled out in view of the low rate constant recently calculated for this reactioiP^. A possibility for 
the origin of water in the innermost region^ of the envelope is pulsationally induced shocks which 
results in a chemical stratification different from thermodynamic equilibrium (TE) chemistry. However, 
IRC +10216 has a C/0 abundance ratio of 1.4IS Recent non-TE calculation^!^ have shown that for 
a carbon-rich star with an even lower C/0 ratio of 1.1, water should be almost completely absent in 
the inner wind, and might only exist between 1 and 1.4 i?^, with a [H2O/H2] peak abundance around 
5 X 10"'^. Simulating this situation and assuming the (too high) water abundance of 5 x 10"^ over the 
full region between 1 and 1.4 i?^, yields water line predictions being a factor 3 to 10 too low compared 
to the PACS and SPIRE observations, ruling out the shock-induced non-TE chemistry as possible 
cause of water. 

An alternative origin for the warm water vapour could be provided by photochemistry in the 
inner regions of the CSE of IRC +10216. For a strictly isotropic and homogeneous mass loss pro- 
cess, the inner regions are well protected against the interstellar ultraviolet (UV) radiation by the 
circumstellar material located outwards (the visual extinction of interstellar light is more than 100 
magnitudes for the innermost regions in IRC +1021(1^. Circumstellar envelopes are, however, not 



perfectly spherical but have inhomogeneities and a more or less clumpy structure. Observational ev- 
idence for the clumpy structure of the envelope around IRC +10216 has been found both at small 
and large scales through observations at near-infrared and visible wavelength^^^] ESI ^s through 

millimeter- wave observations of different molecule^^D EH ES] ^ggg Supplementary Information) . The ex- 
istence of a clumpy structure allows for a deeper penetration of a fraction of interstellar UV photons 
into the inner circumstellar layers, thus promoting dense and warm UV-illuminated inner regions. 
In such an environment water can be formed from the photodissociation of the major oxygen-carrier 
molecules, mostly ^'^CO and SiO (^^CO is hard to photodissociate due to self-shielding effects), and 
the subsequent liberation of atomic oxygen, which then converts into water through the chemical 
reactions: 

O + H2 ^OH + H, and (1) 

OH + H2^H20 + H, (2) 

which are only rapid enough at temperatures above ~300 K. We stress that this mechanism does 
not necessarily extend to the whole inner CSE but only to those inner clumpy regions which are 
more exposed to the interstellar UV field. For the case where 10% of the total circumstellar mass is 
illuminated by interstellar UV photons through a cone where matter only fills 70% of the solid angle 
of arrival of interstellar light, we predict that water would be formed in the inner envelope with a 
maximum abundance relative to H2 in excess of 10~^ (see Fig. 3 and the Supplementary Information). 
The predicted water line strengths for the case of a minor UV-illuminated component are shown in 
green in Fig. 1. The deduced amount of water is 0.003 earth masses. 

AGB stars also have a soft UV stellar radiation fielcP31E51_ Pulsationally induced shocks might 
generate a surplus of UV photons close to the stellar photosphere. However, even for a clumpy inner 
envelope, the densities just above the stellar photosphere are so high that the UV photons will be 
severely attenuated in the first few 10^^ cm, in contrast to the outer envelope, where the material is 
much less dense. 



The penetration of interstellar UV photons will yield the formation of hydrides, other than H2O, 
in the inner envelope through successive hydrogenation reactions of heavy atoms (nitrogen, carbon 
or sulphur). Ammonia (NH3, see Fig. 3) is an interesting example as it has been observed in IRC 
+10216 and also from the CSEs of oxygen-rich AGB stars with abundances relative to H2 in the 
range 10~^— 10~^ (EIIID^ much larger than the 1 x 10~^^ predicted by thermochemical model^^'^. 
Other molecules, such as HC3N, which are typically formed by photochemistry in the outer layers, 
also have increased abundances in the inner regions, as seen in Fig. 3. The predicted higher abundance 
of HC3N in the inner regions is confirmed by our recent observations of HC3N (see Fig. 4) with the 
IRAM 30 m telescope at Pico Veleta, for which the line profiles point towards the existence of a warm 
inner component. The discovery of high-excitation H2O lines in the inner warm and dense envelope 
of an evolved carbon-rich star questions our knowledge of the envelope chemistry and outlines the 
importance of UV induced photochemistry in the CSEs of evolved stars. In the case of oxygen-rich 
environments, the same mechanism predicts high abundances of carbon-rich species, such as HCN, 
CII4 and C^^, as has already been observed in several target^^Sl EQ] 
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Figure 1: Unblended ortho- and para-water lines detected with Herschel in 

IRC +10216. 39 ortho-H20 and 22 para-H20 lines are identified, including low and high- 
excitation lines (see also Supplementary Figs. 1-5). Panel (a) shows the ground-state ortho- 
H2O line as observed with SWA^ (units of intensity, T4, in K versus velocity with respect 
to the local standard of rest, vlsr, being — 26kms~^ for IRC +10216). Panels (b) - (i) show 
the continuum subtracted flux (in Jy) versus wavelength of observation for ortho-H20 lines, 
while panels (1) -(o) show several para-H20 lines. The HCN 1^2 = 4 contribution to the ortho- 
H20(42,3 — 81^2) line in panel (d) is ~15 Jy. 

The coloured lines show the non-local thermodynamic equilibrium predictions (see Supplemen- 
tary Information) for the different chemical mechanisms proposed as cause of water vapour 
in the envelope of IRC +10216. For the red and orange lines, the predictions of ortho-H20 
represent envelope models with a constant fractional abundance of ortho-H20 (relative to H2) 
out to 4 X lO^^cm, where it is photodissociatecP. The abundance of ortho-H20 was derived 
from the SWAS observations (see panel (a)). The red model simulates ortho-H20 originating 
in the intermediate envelope, with inner envelope radius, Rmt, of 2.1 x lO^^cm and the de- 
rived abundance, [ortho-H20/H2], of 2.5 x 10~^. This model applies both for the hypothesis of 
the vaporization of icy bodies and the Fischer- Tropsch catalysis mechanism, predicting water 
around a few times lO^^cm. The orange model assumes ortho-Il20 to be generated in the outer 
envelope, with Rint of 4.3 x lO^^cm and [ortho-H20/H2] equal to 6.7 x 10^^. Finally, the green 
model shows the model predictions for ortho-H20 present in the inner envelope, with a radial 
distribution of the fractional abundance as shown in Fig. 3. The last row shows four para- water 
lines, where the blue line is based on the fractional abundance distribution as show in Fig. 3, 
using an ortho-to-para H2O ratio of 3:1. 



Figure 2: Schematic overview (not to scale) of the envelope around a carbon-rich 
AGB star. Several chemical processes are indicated at the typical temperature and radial 
distance from the star in the envelope where they occur. Shock-induced non-equilibrium chem- 
istry takes place in the inner wind envelop^^, dust-gas and ion-molecule reactions alter the 
abundances in the intermediate wind zone, molecules such as CO and HCN freeze-out at inter- 
mediate radiP^, and the penetration of cosmic rays and ultraviolet (UV) photons dissociates 
the molecules and initiates an active photochemistry that creates radicals in the outer wind 
region!^. The different mechanisms hitherto proposed as origin of water in a carbon-rich envi- 
ronment are indicated at the bottom of the figure in grey at the typical distances where they 
occur. The penetration of the interstellar ultraviolet photons in a clumpy circumstellar envi- 
ronment is shown in blue. The resulting chemical processes important for the creation of H2O 
and HC3N are indicated in red. 



Figure 3: Fractional abundances in the clumpy circumstellar environment of 

IRC +10216. The chemical model simulates a clumpy envelope structure, where a fraction 
of the interstellar UV photons is able to penetrate deep into the envelope (see Supplementary 
Information). This figure shows the predictions for the radial distribution of the fractional 
abundances relative to H2 of H2O, NH3, and HC3N for a model with a minor UV-illuminated 
component superposed on a major UV-shielded component. The minor component shown in 
this figure contains 10% of the total circumstellar mass {Jm = 0.1), which is illuminated by 
interstellar UV photons through a cone where matter fills 70% of the solid angle of arrival of 
interstellar light {fn = 0.3). Dashed lines correspond to the minor UV-illuminated component, 
continuous lines to the major UV-shielded component, and thick grey lines to the weighted 
average abundance over the two components. The weighted average abundance is computed as 
Tiir) = (1 - /M)Xr^^°'(r) + {fM)X~{r), where X°'^j°''(r) and X~{r) are the abundance 
of the species i in the major UV-shielded and minor UV-illuminated component, respectively, as 
a function of radius r. Note that for H2O and NH3 the contribution from the major UV-shielded 
component is negligible, and the dashed and thick grey line coincide. 



Figure 4: HC3N as observed in the envelope of IRC +10216. The HC3N J=29 28, J=33 
32 and J=37-36 hnes have been observed with the IRAM telescope. All three HC3N lines show 
a clear flat-topped profile. The HC3N J=29-28 line is shown in this figure, and its line profile 
is compared to that of the AlF J=8-7 line (black line). Both lines were observed with the same 
telescope, the same beam, the same pointing, and have been calibrated in the same way. While 
the AlF line profile is U-shaped, the HC3N line is fiat topped, clearly indicating that HC3N 
arises from gas extending to inner radii. The grey lines show two model predictions: the dashed 
grey line corresponds to a model prediction only taking the major UV-shielded component into 
account, the full gray line shows the theoretical line profile for a model including both the major 
UV-shielded component and the minor UV-illuminated component. The spectrum is plotted 
in terms of intensity in K) versus velocity with respect to the local standard of rest (vlsr 
in kms~^). The LSR stellar radial velocity of IRC +10216 is — 26kms~^. 
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1. Radiative transfer 

PACS and SPIRE observations of IRC +10216 have been obtained in November 2009^3. The 
PACS data (55 - 200 yum) presented in the first paper describing the Herschel observations^'^ 
were calibrated using a ground-based absolute flux calibration. Recently, a new absolute flux 
calibration for the PACS data has been obtained^^, yielding a decrease in absolute line flux of 
30% in the blue wavelength range (A < 98 yum) and of 10% in the red wavelength range. The 
estimated calibration uncertainty on the PACS line fluxes is ~30%, but is higher for wavelengths 
shortward of 63 fim due to non-optimal detector flatfielding. The PACS wavelength accuracy 
is around 20%^^. For the SPIRE data (194 - 670 /im) the absolute flux calibration uncertainty 
is 15-20% between 194 and 313 /im, and 20-30% between 303 and 500 /im, going up to 50% 
beyond 500 /im. 

These newly calibrated data (see Supplementary Figs. 1-5) have been used to rederive 
the thermal structure of the envelope around IRC +10216 through non-local thermodynamic 
equilibrium (non-LTE) radiative transfer modelling of the ^^CO lines, using the GASTRoNOoM 
code^^EH. The PACS and SPIRE ^^CO lines cover energy levels from J = 3 (at 31 K) to J = 47 
(at 5853 K) and trace the envelope for radii r ^ 1 x 10^^ cm (r < 2000 i?*). The 60 lowest 
levels in the ground and flrst-vibrationally excited state are included in the radiative transfer 
modelling^^. The velocity structure was obtained by solving the momentum equation. The 
(circum) stellar parameters are given in Table 1. Speciflcally, we derive a gas mass- loss rate of 
2.1 X 10"'^MQ/yr for a [CO/H2]-ratio of 1 x 10~^*E1) ^n^j ^ kinetic temperature that is given 
by r(r) = Tcff(-R^) x (r/i?^)~°'^^ (for r < 10 i?^) and T(r) oc r""-^ beyond that radius. The 
estimated uncertainty on the mass-loss rate is a factor of 2. The derived kinetic temperature 
agrees very well with the results obtained for the inner envelope from a detailed modelling of 
C2H2 and HC^P. 



Using the derived thermodynamical structure, the level populations for ortho- and para- 
H2O lines are calculated using the GASTRoNOoM-code for different stratiflcations of the water 
abundance throughout the envelope. The radiative transfer modelling includes the 45 lowest 
levels in the ground state and flrst vibrational state (i.e. the bending mode 1^2 = 1 at 6.3 /im). 
Level energies, frequencies and Einstein A coefficients are extracted from the Barber water line 
lislP^. The H2O-H2 coUisional rates are taken from the Faure databas^^. The effect of including 
excitation to the flrst excited vibrational state of the asymmetric stretching mode {u^ = 1) has 
been tested, and is found to be negligible. 
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Supplementary Figure 1: Herschel spectrum of IRC +10216 compared to water line 
predictions. The Supplementary Figures 1-5 show the continuum subtracted Herschel PACS 
and SPIRE spectrum of IRC +10216 (in black). The coloured lines are the model predictions 
for ortho- and para-water with parameters as described in the caption of Fig. 1 in the main 
journal. The detected water lines and the ^^CO and ^^CO lines are indicated. In total, 39 
ortho- water lines and 22 para-water lines are detected in the PACS and SPIRE spectrum. Most 
of the water lines seen in the Supplementary Figures 1-5 are blended with another molecular 
emission line arising from HCN, SiO, SiS, CO, etc.; a sub-sample of (almost) unblended water 
lines is presented in Fig. 1 in the main journal. 
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Supplementary Figure 2: see caption of the Supplementary Figure 1. 
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Supplementary Table 1: (Circum) stellar parameters of IRC +10216. T^s is the effective 
stellar temperature, i?* denotes the stellar radius, Voo the terminal velocity of the wind, M the 
gas mass-loss rate, and -Rdust the dust condensation radius. 



Teff[K] 




2330® 


M[Mo/yr] 


2.1 X 10"^ 


R. [10^3 


cm] 


5.1<E1I) 


-Rdust [R-k] 


5.6(113) 


[CO/H2] 


[10-3] 






30(13) 


distance 


[pc] 


180® 


Voo [kms-^] 


14.5(E2I) 



Continuum emission provides the dominant excitation source for H2O, mainly through 
absorption of 6 /xm photons in the 1^2 band and subsequent decay^. The spectral energy dis- 
tribution (SED) from 2 fim to 1 mm has been modelled using the MCMAX cod^^, a radiative 
transfer code specialised for the treatment of optically thick dusty circumstellar media. The rel- 
evant stellar and circumstellar input parameters are identical to the ones used in the molecular 
emission line modelling. The best fit to the ISO SWS and LWS (Infrared Space Observatory, 
Short Wavelength Spectrometer and Long Wavelength Spectrometer) spectrophotometric data 
between 2.4 and 197 /im, supplemented with photometric data point^, requires the inclusion 
of amorphous carbon, iron, silicon carbide and magnesium sulfide as dust species. Assuming 
spherical grains, their respective derived mass fraction abundances are 0.72, 0.05, 0.13 and 0.10. 
From these relative abundances, the mean specific density of the dust grains, the dust mass-loss 
rate and the wavelength-dependent absorption efficiencies are calculated to provide the dust 
component during the modelling of the molecular emission lines. The derived dust mass-loss 
rate is 8 x 10~^ Mq/jt, yielding a dust-to-gas ratio of ~ 4 x 10-^. Estimated uncertainties are 
a factor of 2. The mean specific dust density is 2.41 gcm^^^ 

The water fractional abundance as derived from the SWAS observations (see Fig. 2) is 
slightly higher than results presented already in the literature^. This mainly reflects a difference 
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Supplementary Figure 3: see caption of the Supplementary Figure 1. 
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Supplementary Figure 5: see caption of the Supplementary Figure 1. 
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Supplementary Figure 6: see caption of the Supplementary Figure 1. 



in (circum)stellar parameters and dust opacities. Using the same input parameters, similar 
abundances are derived. 



2. Chemical modelling 



There is a lot of evidence for a clumpy structure of the envelope around IRC +10216 from 
astronomical observations carried out at small and large scales in the visible and near-infrared 
wavelength ranges^ These optical and near-infrared observations trace the dust density 
distribution, and thus the gas density distribution if both dust and gas are position coupled^". 
Millimetre- wave interferometric observations have mapped the spatial distribution of CO, find- 
ing a large scale clumpy structure^^. Other molecules, such as CN, C2H, C4H, and HC5N show 
a clumpy and hollow quasi-spherical distribution with a lack of molecular emission in the NNE 
to SSW direction^"'^^. The hollow distribution can be interpreted in terms of photochemistry 
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of the parent species, creating reactive species which then undergo further chemical reaction^. 
The conical holes in the NNE and SSW directions have been interpreted in terms of a molecular 
density minimum in that part of the envelope^^, or the fact that IRC +10216 is just leaving the 
AGB and that a bipolar outflow has dug out two symmetric conical holes at a position angle 
PA~20°'S2I). The existence of these clumps and conical holes opens a path for interstellar UV 
photons to penetrate deep into the envelope. 

Whilst the existence of clumpy structures and cavities in the envelope of IRC +10216 is 
not in doubt according to the observations, it is difficult to model this complicated structure, 
since the different parameters characterizing the three-dimensional clumpy structure are not 
well constrained by the observations. We have thus adopted a simple modelling approach to 
simulate the clumpy structure and to investigate the effects on the circumstellar chemistry. 
We considered that the envelope consists of two different components: a major one which is 
well shielded against interstellar UV photons and only starts to undergo photochemistry in the 
outer layers, and a minor one (which accounts for a fraction J'm of the total circumstellar mass) 
for which the shielding material located in the radial outward direction is grouped into clumpy 
structures leaving a fraction J'q of the solid angle of arrival of interstellar light free of matter. 
We name these two components 'the major UV-shielded' and 'the minor UV-illuminated' com- 
ponent. As might be clear the minor UV-illuminated component is not completely exposed 
to UV photons, but is also partially shielded against UV photons by the clumps located on 
the way to the interstellar medium. Likewise, the major UV-shielded component starts to be 
partially exposed to UV photons in the outer envelope. 

The physical conditions in both components are assumed to be the same, except that for 
the minor UV-illuminated component the visual extinction against interstellar light is reduced, 
as compared to that of the major UV-shielded component, due to the empty space between 
the clumpy material. The observations quoted in the previous paragraph have not provided 
an estimate of the density in the interclump medium. In our chemical model we have assumed 
that the density in the interclump medium is formally zero (in practice it is low enough to allow 
interstellar UV photons to penetrate without significant extinction). The dumpiness is assumed 
to extend down to the innermost regions of the envelope (r ~10^^cm). Therefore, besides the 
standard physical parameters of any model describing a spherical circumstellar envelope (see 
Supplementary Table 1), we utilize two further phenomenological parameters to take care of 
the effects of dumpiness on the chemistry: for the minor UV-illuminated component J'm is the 
fraction of the total circumstellar mass in the minor component and fu is the fraction of the 
solid angle of arrival of interstellar light which is free of matter. 

The gas phase chemical composition of both the major UV-shielded and the minor UV- 
illuminated components are obtained from chemical kinetic calculations starting at the inner- 
most layer (r = 10^^ cm) up to a distance of r = 10^^ cm, where most molecules have already 
been destroyed by interstellar UV photons. The chemical network utilized is essentially that 
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used in previous work on the chemistry of oxygen in IRC +10216^, with updates for some reac- 
tion rate constants and photodissociation/ionization rates from the recent hteratur^^l^. The 
interstellar UV field has been decreased by a factor of 2 with respect to the standard interstellar 
UV fielcP^, to predict the angular position of the maximum abundance of C2H and CN at 15" 
and 20" from the star, respectively, in agreement with interferometric observation^. We note 
that when fitting the peak abundance distribution of these molecules, there exists a degeneracy 
between various physical parameters, such as the interstellar UV field, the mass loss rate, the 
extinction law, and the distance to IRC +10216. Furthermore, the exact value of the local 
interstellar UV field around IRC +10216 is not precisely known. Recently, an extended UV 
structure surrounding IRC +10216 has been discovered at a radial distance of about 12'®. 
This structure has been attributed to the interaction of the stellar wind with the interstellar 
medium, and it might contribute to the ambient UV field seen by the circumstellar envelope 
of IRC +10216. The exact value of the interstellar UV field around IRC +10216 is, however, 
not critical in determining the flux of UV photons reaching the inner layers of the minor UV- 
illuminated component. What really counts is the degree of dumpiness, i.e. the parameter 
Jq. a more detailed description of the chemical modelling approach for this kind of environ- 
ment, with a major UV-shielded and a minor UV-illuminated component, for both carbon- and 
oxygen-rich circumstellar envelopes with different mass-loss rates will be presented in another 
article^^. 

For the envelope around IRC +10216, the inner regions of the major UV-shielded com- 
ponent (surrounded by a smooth envelope) are indeed well shielded against interstellar UV 
photons, with a visual extinction larger than 50 mag for the regions down to 5 R^. Assuming 
that the minor UV-illuminated component is illuminated by interstellar UV photons in a cone 
where matter is just filling 50-80% of the solid angle of arrival of interstellar light (/^ = 0.5- 
0.2), results in a substantial decrease of the effective visual extinction in the inner envelope from 
more than 50 mag to less than 0.5 mag, which indeed allows for photochemistry to take place 
in the inner envelope. In the final model, whose results are shown in Fig. 3 of the main article, 
we have adopted a value of 0.3 for Jq and assumed that the minor UV illuminated component 
accounts for 10% of the total circumstellar mass (/a/ = 0.1). The resulting abundance radial 
profile for each species i are then obtained as the weighted-average of the abundance in both 
components, i.e. as X^r) = (1 - fM)X^''^°\r) + (/jv/)X™(r), where X^^^°\r) and X™(r) 
are the abundance of the species i in the major UV-shielded and minor UV-illuminated compo- 
nent, respectively, as a function of radius r. This approach assumes that the major UV-shielded 
and the minor UV-illuminated components are not distributed over a given preferred radial di- 
rection but over many different radial directions, in such a way that we can average over all 
directions to get a mean value of the abundance Xi at each radius r. Whilst it is clear that 
in reality both components cannot coexist in the same radial directions, this approach permits 
us to describe the distribution of the molecular abundances under spherical geometry, which 
facilitates the radiative transfer modeling described in Sect. 1. We note that the effect of a 
clumpy medium on spectral line formation can be characterized by a single number multiplying 
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the mean column opacitj^. Even for an arbitrarily complex clump distribution, the effect can 
be reduced to a single scaling number. The scaling factor reduces the effective opacity of a 
clumpy medium compared to a smooth gas of the same mean column density. The scaling 
factor is directly retrieved from the chemical modelling (see Fig. 3 in the main article) and 
taken into account when modelling the H2O line profiles. 

As shown in Fig. 3 of the main article, water vapour reaches a weighted-averaged abun- 
dance relative to H2 in excess of 10~^. The photodissociation radius for each molecule is 
calculated from the chemical network taking the two components into account. According to 
the model water extends from the innermost regions (r ~ 10^^ cm) up to a radial distance of 
(1-3) X 10^^ cm. In this region water vapour keeps a roughly constant abundance due to the 
balance between the formation processes (the reaction of O and H2 followed by that of OH and 
H2) and the destruction processes (photodissociation by the ambient interstellar UV field). The 
photodissociation radius of water, (1-3) x 10^^ cm, is noticeably lower than that predicted by 
other hypotheses, 4 x 10^^ cm (see main article). The smaller photodissociation radius obtained 
by our model is fully consistent with the SWAS, PACS, and SPIRE observations, as shown in 
Fig. 2 of the main article. 
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